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Abstract 

We review the theoretical approaches for investigations of the high energy hadron-hadron 
scattering in the Regge kinematics. It is demonstrated, that the gluon in QCD is reggeized 
and the Pomeron is a composite state of the reggeized gluons. Remarkable properties of the 
BFKL equation for the Pomeron wave function in QCD and supersymmetric gauge theories 
are outlined. Due to the AdS/CFT correspondence the BFKL Pomeron is equivalent to 
the reggeized graviton in the extended N=4 SUSY. The properties of the maximal tran- 
scendentality and integrability are realized in this model. The BDS multi-gluon scattering 
amplitudes are investigated in the Regge limit. They do not contain the Mandelstam cuts 
and are not valid beyond one loop. It is shown, that the hamiltonian for these composite 
states coincides with the hamiltonian of an integrable open Heisenberg spin chain. 

1 High energy interactions 

Hadron-hadron scattering in the Regge kinematics 

s = (pa + PB? = (2E) 2 »f = -{p A , - p A ) 2 ~ m 2 (1) 

is usually described in terms of a i-channel exchange of the Reggeon (see Fig.l) 

• m , e -* v M*) + v 

A p (s, t) = £ p (t) g(t) 8*® g(t) , j p (t) = j + a't , £ p (t) = - F , (2) 

where j p (t) is the Regge trajectory which is assumed to be linear, jo and a' are its itercept and 
slope, respectively. The signature factor £ p is a complex quantity depending on the Reggeon 
signature p = ±1. 



A special Reggeon - Pomeron with vacuum quantum numbers is introduced to explain an 
approximately constant behavior of total cross-sections at high energies and a fullnllment of 
the Pomeranchuck theorem cr^/ahh — * 1- Its signature p is positive and its intercept is close 
to unity j£ = 1 + A , A < 1. 

The particle production at high energies can be investigated in the multi-Regge kinematics 
(see Fig. 2) 
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s > si , s 2 , ... , S n+ i » ti , t 2 , ••• , tn+1 , (3) 

where s r are squares of the sum of neibouring particle momenta k r -±, k r and —t r are squares 
of the momentum transfers q r . The production amplitude in the framework of the Regge model 
is also expressed in terms of the Reggeon exchanges in each of t r -channels 

n+1 

A^ 2+n ~ n 4 p( * r) , (4) 

r=l 

where we neglected the signature factors, which will be discussed later. 



2 Gluon reggeization in QCD 

In the Born approximation of QCD the scattering amplitude for two colored particle scattering 
is factorized 

M AB { s i t)\ Bo rn = ^A'A ~T ^B'B > ^A'A = 9 T A'A <^A A ,A A , (5) 

where T° are the generators of the color group SU(N C ) in the corresponding representation and 
A r are helicities of the colliding and final state particles. In the leading logarithmic approxima- 
tion (LLA) the scattering amplitude in QCD can be written as follows pQ 

M% B B ' (s, t) = MH' (s, t)\ Born s"® , a s In s ~ 1 , (6) 

where the gluon Regge trajectory is 

2 a s N c 2 f tPtf-^k _ ( \q 2 \ 1\ 

Here the extra dimensions 2e — > —0 were introduced to regularize the infrared divergency, the 
parameter fi is the renormalization point and we used the notation 

2vr v ; ' 

where 7 is the Euler constant 7 = — "0(1). This Regge trajectory was calculated also in two-loop 
approximation in QCD [2] and in supersymmetric gauge theories [3]. 

Further, the gluon production amplitude in the multi-Regge kinematics at LLA can be 
written in the factorized form [1] (see Fig. 2) 

g^l g^2 s LO n 

M 2 ^ 1+n = 2 S T% AT ^gT^ Cl C(q 2 , qi ) C(q n , g n _i) T% B . (8) 

mr ml mr 

The Reggeon-Reggeon-gluon vertex for the produced gluon with a definite helicity is 

C(q2,qi) = -^j, (9) 

q 2 ii 
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where we used the complex notations for the transverse components of particle momenta. 
It gives a possibility to calculate the total cross-section [1] 



[ dT n \M 2 _ 1+n \ 2 , (10) 
where T n is the phase space for the produced particle momenta in the multi-Regge kinematics. 

3 BFKL equation 

Because the production amplitudes in QCD are factorized, one can write a Bethe-Salpeter-type 
equation for the total cross-section at- Using also the optical theorem it can be presented as the 
equation of Balitsky, Fadin, Kuraev and Lipatov (BFKL) for the Pomeron wave function [lj 

EV{p u fc) = H 12 p 2 ) , A = -^E, (11) 



where at ~ Smax an d the BFKL Hamiltonian in the coordinate representation p is 

#12 = In \pi P 2\ 2 + -^(hi \pn\ 2 )piP* 2 + -^-(h \pi2\ 2 )plP2 - 4^(1) (12) 

PlP * p\p 2 

and p\2 = pi — P2- It is invariant under the Mobius transformations [HE] 

apk + b 

Pk -» —7 13 

cpk + d 

and has the property of the holomorphic separability 

Hu = h u + h* 2 , hi 2 = ln(pip 2 ) + — ^(pu) Pi + — ln(pi 2 ) P2 - 2-0(1) • (14) 

Pi P2 

Here we used the complex notations for two-dimensional transverse coordinates and their canon- 
ically conjugated momenta. The conformal weights for the principal series of unitary represen- 
tations of the Mobius group are 

m = 7 + n/2 , m = j — n/2, r y = l/2 + iv, (15) 

where 7 is the anomalous dimension of the twist-2 operators and n is conformal spin. 

The Bartels-Kwiecinski-Praszalowicz (BKP) equation for colorless composite states of sev- 
eral reggeized gluons has the following form [6] 

E^>(p 1 ,...)=H^(p 1 ,...) , H = Y,™H kl , (16) 

k<l c 

where H^i is the BFKL hamiltonian. Apart from the Mobius invariance its wave function in 
the multi-color QCD (N c — > 00) has the property of the holomorphic factorization [7] 

*(pi,...,p n ) =J2 a r,s^r{pi,-,Pn)^s(pl,-,Pn), (17) 
r,s 

where the sum is performed over a degenerate set of solutions for the corresponding holomorphic 
and anti-holomorphic equations. The BKP equation has the duality symmetry p k — > pk,k+i —* 
Pk+i (k = 1, 2, ...,n) [8] and n integrals of motion q r , q* 0. The corresponding hamiltonians 
h and h* are local hamiltonians of the integrable Heisenberg spin model, in which spins are 
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generators of the Mobiuos group [10J. We can introduce the transfer (T) and monodromy (t) 
matrices according to the definitions [9 J 

n 

T(u) =trt{u) = ^u n - r q r , t(u) = L x L 2 ...L n , (18) 

r=0 

L k =( U + P 2 kPk Pk ), t{u) = ( j^f) B J-;\ ) . (19) 
V -PkPk u-pkPkJ' KJ \C(u) D(u) J v ' 

The matrix elements of t(u) satisfy some bilinear commutation relations following from the 
Yang-Baxter equation [9] 

tp (u) tp (v) (v-u) = l s J s J (v - u) 4 (?) 4 (u) , l(u) =ul + iP, (20) 

'1 '2 *1*2 

where l(u) is the monodromy matrix for the usual Heisenberg spin model and P is the per- 
mutation operator. This equation can be solved with the use of the Bethe ansatz and the 
Baxter-Sklyanin approach I12j. 



4 Pomeron in TV = 4 SUSY 

One can calculate the integral kernel for the BFKL equation also in two loops [13] . Its eigenvalue 
can be written as follows 

u = 4 a x (n, 7) + 4 a 2 A(n, 7) , a = g 2 N c /{\^ 2 ) , (21) 

where 

X (n, 7) = 2* (1) - * (7 + \n\/2) - * (1 - 7 + \n\/2) (22) 

and \F(x) = T'{x)/T{x). The one-loop correction A(n,7) in QCD contains the non-analytic 
terms - the Kroniker symbols 5| n |,o an d <V|,2 P]- But in N = A SUSY they are cancelled and 
we obtain for A (71,7) the following result in the hermitially separable form [3l 114] 

A ( „, 7) = m) + «m - m±f<m , M _ 7 + m , (23) 



p{M)=P>(M) + \cm^\z) = \ 



(24) 



It is interesting, that all functions entering in these expressions have the property of the maximal 
transcendentality [13] • In particular, (f>(M) can be written in the form 

4>{M) = 3C(3) + " (M) - 2$(M) + 2/3' (M) - *(M)) , (25) 

$(M ) = f; |4- (*< ( * + 1) - SM-ffl) . (20) 

Here ^(M) has the transcedentality equal to 1, its derivatives have transcedentalities n+1, 
the additional poles in the sum over k increase the transcedentality of the function <J>(M) up to 
3 being also the transcendentality of C(3)- The maximal transcendentality hypothesis is valid 
also for the anomalous dimensions of twist-2 -operators in N = 4 SUSY |15[ [16] contrary to the 
case of QCD [IT]. 

The stationary BFKL equation in the diffusion approximation can be written as follows [1] 

j = 2-A-Dv 2 , (27) 
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where v is related to the anomalous dimension 7 of the twist- 2 operators [13] 

1 = l + l^ + iu. (28) 

The parameters A and D are functions of the coupling constant a and are known up to two loops. 
Higher order perturbative corrections can be obtained with the use of the effective action [181119]. 
For large coupling constants one can expect, that the leading Pomeron singularity in N = 
4 SUSY is moved to the point j = 2 and asymptotically the Pomeron coincides with the 
graviton Regge pole. This assumption is related to the AdS/CFT correspondence, formulated 
in the framework of the Maldacena hypothesis claiming, that N = 4 SUSY is equivalent to the 
superstring model living on the 10-dimensional anti-de-Sitter space \20 \ \21 \ [22] . Therefore it is 
natural to impose on the BFKL equation in the diffusion approximation the physical condition, 
that for the conserved energy-momentum tensor 0^ LU {x) having j = 2 the anomalous dimension 
7 is zero. As a result, we obtain, that the parameters A and D coincide [16J. In this case one 
can solve the above BFKL equation for 7 

J = U-2)( 1 - - 1/A V (29) 

} \2 l + Vl + (i-2)/Ay 

Using the dictionary developed in the framework of the AdS/CFT correspondence [21], one 
can rewrite the eigenvalue relation for the BFKL kernel in the form of the graviton Regge 
trajectory [16] 

j = 2 + yt, t = E 2 /R\ <*' = Y A - ( 3 °) 

On the other hand, Gubser, Klebanov and Polyakov predicted the following asymptotics of the 
anomalous dimension at large a and j [23] 

1/4 



7|a,i^oo = = V^j a 1/4 • (31) 

As a result, one can obtain the explicit expression for the Pomeron intercept at large coupling 
constants [161121] 

j=2-A, A = i-a" 1 / 2 . (32) 

Note, that in Ref. [25] it was argued, that for N = 4 SUSY the evolution equations for 
anomalous dimensions of quasi-partonic operators are integrable in LLA. Later such integra- 
bility was generalized to other operators [2B] and to higher loops [27]. Using additionally the 
maximal transcendentality hypothesis the integral equation for the so-called casp anomalous 
dimension was constructed in all orders of perturbation theory [281 129| . Further, the anoma- 
lous dimension of twist- 2 operators in four loops was calculated [30J, but due to the absence 
of so-called wrapping contributions in the asymptotic Bethe anzatz the obtained results do not 
agree with the BFKL predictions (3J E] . 



5 Bern-Dixon-Snirnov scattering amplitudes in N = 4 SUSY 

To calculate higher order corrections to the BFKL equation in QCD and supersymmetric models 
one should know production amplitudes in higher orders of perturbation theory. Several years 
ago Bern, Dixon and Smirnov suggested a simple anzatz for the multi-gluon scattering amplitude 
with the maximal helicity violation in the planar limit aN c ~ 1 for the N = 4 super-symmetric 
gauge theory [31J. It turns out, that this amplitude is proportional to its Born expression. The 
proportionality coefficient M n for n external particles is a function of relativistic invariants and 
can be written as follows 

00 

In M n = a 1 (/« (e) (/« (k) + F« (0)) + C« + (e)) , (33) 
1=1 
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oo oo oo 

/%) = /<*> + efP + e 2 /f , 7(«) = 4 £ a'jf , /3(a) = £ a 1 /? , 8(a) = £ a'/? , (34) 

Z=l i=l !=1 

where En\e) can be neglected for e — > and the constants and fr\e) are known up to 
rather high order of the perturbation theory. In particular, 7(a) is the so called cusp anomalous 
dimension which was found in all orders [281 [29] 

7 (a) =4a-4C 2 a 2 + 22C 4 a 3 + .... (35) 
The singular function In\e) is given below 



S"M = -5?E(: 



M 2 



" s i,i+l 



(36) 



and the finite remainders Fn are expressed in terms of logarithms and dilogarithms. 

In Ref. [32] the BDS anzatz was investigated in the Regge kinematics (see also Ref. |33j). 
In particular, the elastic amplitude has the Regge asymptotics 

M2-.2 = T(t) (^y* T(t) (l + <5(e)), (37) 

where /i 2 is the renormalization point, 

7(a) , -t f a da' / 7 (a') , , 
4 ^ 2 J a' V 4e ' 



C 2 (ln.. 2 



a 2 + ... (38) 



^ 2 2ey 2 

is the all-order gluon Regge trajectory obtained from the BDS formula |32| and 



, -* 1 , 
In — k + - I a + 



lnr(t) = In- / _(^-Z + ^) + ^ + ^_Z C2 







da' 


/^ 2 J 





~V 


a da' 


In 


a 








a' 



(39) 



is the vertex for the Reggeized gluon coupling to the external particles. Note that the pertur- 
bative expansion for u{t) is in an agreement with its direct calculations done initially in the 
MS-scheme g|. 

One can verify that in all physical regions the BDS amplitude for one gluon production in 
the multi-Regge kinematics can be obtained with the use of an analytic continuation from the 
expression [32] 

M 2 ^ 3 /_, \ *W*»)/_ ()! \«(0) /_ S2 Wfc)-«(*0 /-sk\ U{tl) , . 

ci + — »- —3- c 2 , (40) 



r(*i)r(t 2 ) w y W y V /^ 2 7 V ^ , 

where the coefficients c, are real 

i„/ , m sin it ( uj ( t{) — d>r) . 
= r(t 2 ,t!,ln-«) . V ; 1; ; r; 41 
sm7r(w(iij — w(t 2 Jj 

= rfatx.ln-zc) \; 2; J r; 42 

sm7r(w(i 2 J — w(£iJJ 

Here (fir is the phase of the Reggeon-Reggeon-gluon vertex T, i.e. 

r(t2,ti,ln«-t7r) = |r(t 2 ,*i,ln-«)| e i7r< ^ , (43) 



1. 



defined by the expression 

, „. , 7(a) , 2 -k 1 f a da' a (Ma!) f3(a') e . ,,\ 

lnr (t 2 ,hM-K) = -^ln 2 — --/ —]n-lJ±-J. + 2U. + 6{a')) 

16 n A 2 ,/ a' a' \ 4e z e J 

In the above dispersion-type representation for the production amplitude we can use the 
reality condition for the produced gluon 

K _> flf 2 . = k\ , (45) 
s 

where k± is the transverse component of its momentum (k±pA = k±ps = 0). 

In a similar way two gluon production amplitude in the multi-Regge kinematics almost in all 
physical regions can be obtained by an analytic continuation from the following dispersion-like 
representation for the BDS expression 



Tir / \ U)(tl)—U>(t2) / \ Ul(t2)— w(<3) / \ LoUs) 
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\ Ul(t2)—U>(tl) / \ w(<3) / \ U){tz) 

\ U)(<2)— w(<3) / \ tofts) — d)(ti) / \ Ul(tl) 



+ ( _^e„-»<« (™y<« 4 (46) 



a* 2 7 V /" 2 



with the real coefficients dj=i,2,3,4,5 satisfying the relations 

di = ci(t 2 ,ti, K12) ci((t 3 ,t 2 , k 23 ) , 
d-2 = C 2 ((t2,h, K 12 ) c 2 ((t 3 , t 2 , K23) , 
d 3 + d 4 = c 2 ((t2,£i, K12) ci((t 3 ,i 2 , k 23 ) , 

^5 = ci((i 2 ,ii,Ki 2 )c 2 ((i 3 ,t 2 ,K 23 ) , (47) 

where 

«12 = (91 - 92)± , k 23 = (<f 2 - g 3 )l • 

However, in the physical kinematical region, where s,s 2 > but Si,s 3 < the Regge 
factorization for the BDS amplitude is broken 

M 2 _> 4 



r(h)T(t 



3) 

-siT** , , . , f-s 2 Y it2) , , , .Js^ {tz) 



C \~^) r(t 2 ,*i,lnKi2-i7r)^-^J r(t 3 ,t2,lnK 23 -ivr) (^-^ ) (48) 



where the coefficient C is given below 

C = exp 



(h + fc 2 )V 



2„2 e 



(49) 
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Similarly for the BDS amplitude describing the transition 3 — > 3 in the physical region, where 
s, s 2 = t 2 > but si, S3 < we obtain the result 



M 3 _ 



r(ti)r(t 3 ) 



c" 



/^ 2 



tj(ti) 



r(t 2 ,ii,lnKi2 + iir) 



-52 



r(t 2 ,ti,ln«23 + in) 



S3 
M 2 



W(*8) 



where the phase factor C" is 

C" = exp 



7^(a) / . n , (9l - <?2) 2 (<?2 - (ft) 2 
( — Z7TJ In — „ 

4 (9i + 93-92) 2 ^ 



(50) 



(51) 



which also contradicts the Regge factorization. The reason for these drawbacks is that just in 
these kinematical regions the amplitudes A 2 ^4 and ^3^3 should contain the Mandelstam cuts 
in the j-pane of the i2-chanel [32] . Therefore the BDS amplitudes for these processes are not 
correct beyond 1 loop. 



6 Mandelstam cuts in the adjoint representation at LLA 

The Mandelstam cuts in the elastic amplitude appear only in the non-planar diagrams because 
the integrals for the Sudakov variables a = 2kP^/s and (5 = 2kps of the reggeon momenta k 
and q — k should have the singularities above and below the corresponding integration contours. 
For the case of the planar diagrams this Mandelstam condition is fulfilled only for inelastic 
amplitudes starting from six external particles in the kinematical region where s,s 2 > and 
si, s 3 < 0. Two reggeons in the i2-channel with an adjoint representation of the gauge group 
SU(N C ) can also scatter each from another. The corresponding contribution to the imaginary 
part in the S2-channel for the amplitude A 2 ^ can be written as follows [32] 

^ m r +ico duj f s2 Y j( \ 

-% 2 M 2 ^ = s 2 J - ) M») (52) 



where the reduced partial wave f 2 {oj) is given by 

h{u) = & € \q 2 \ 2 J cP~ 2e kcp- 2e k' ^ 1 (k,q 2 ,q 1 )G ul (k,k l ,q 2 )^ 3 (k l ,q 2 ,q 3 ) . (53) 
Here $12 are impact factors 

The Green's function G u (k,k' ,q 2 ) satisfies the BFKL-type equation 

^>(*.v.») - (2 ^' ( f + y + vwf^ ( K 8 Gg ") <*■ k '- « ] ' <55) 

where 



K^\k,k';q 2 ) 

5^(k - k') n-\k\ 2 ) + .(-|,2 - k\ 2 ) - 2u(-\q\ 2 )) + ^-k)k\q 2 - kr + ,c. 



(56) 
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The infrared divergencies are extracted in the form of the Regge factor ^ and coincide 
with those of the BDS amplitude, as it should be. The partial wave f 2 (oj) contains the infrared 
divergency only in the one loop 

- I |2 f j2-2e h 1 k Q3 _ a (, I 1 2 I g3 1 2 1\ 

" £ 1921 J q*(k + A*)* \k\z\q 2 - k\* ga(jfe - fc 2 ) ~ 2 V |h + MV J ' 1 j 

which is also compatible with the BDS result. But in the upper loops the iteration of the above 
equation leads to terms which are absent in the BDS amplitude. For example, in two loops we 
obtain for the imaginary part of A 2 ^4 in the S2-channel the following expression 

A ° 2 , , \(H-Qz\ 2 \i2? , \qi - qz?\q2? ,- a . 

A s = — In s 2 In — : — — In — : — — . (58) 

2 \qi\ 2 \k 2 \ 2 \q3r\h\ 2 

It is symmetric with respect to the simultaneous transmutation 

ki k 2 , q\ <-> -q3 ■ (59) 

The same expression is valid for the imaginary part in the s-channel. 

In a similar way we can calculate the s-channel imaginary part of the amplitude for the 
transition 3 — ► 3 

,3^3 a 2 1 j , Igg ~q\ — «?3 1 2 1<?2 1 2 . |g2 — «?i — g3| 2 |<?2| 2 lR „, 

A° s = — In t 2 In 1, l9l , l9 In ; — ^ — ^ . (60) 

2 |fci| 2 |/c 2 | 2 |?3| 2 kl| 2 

Moreover, the BFKL equation for the state with adjoint quantum numbers can be solved ex- 
plicitely and we obtain for the imarginary part in S2-channel [34J 

"~I^L^{wJ ImTJ exp M ,, n) ln S2 ), (01) 

where the eigenvalue of the reduced BFKL kernel for the adjoint representation is 

(\n\ 71 \ 

i>{™ + + W-w + 2V ' (1) ) ' (62) 

It turns out, that the leading singularity of the t2-partial wave corresponds to n = 1 and is 
situated at 

j — 1 = u>(t) + a(41n2 — 2). 



7 Multi-reggeon Mandelstam cuts 

Let us consider now the Mandelstam cuts constructed from several reggeons [35]. The non- 
vanishing contribution from the exchange of r + 1 reggeons appears in the planar diagrams only 
if the number of the external lines is n > 2r + 4. For the inelastic transition 2 — ► 2 + 2r with 
the initial and final momenta pa,Pb and pa', ki,k 2 ,..., k 2r ,pB>, respectively, (see Fig 2) the cut 
exists in the crossing channel with the momentum 

r 2r r+1 

q=PA~PA' -^2k=PB> ~PB+ ^2 kl = ^2 q l ' ( 63 ) 
1=1 l=r+l 1=1 

where q[ are momenta of r + 1-reggeons. The corresponding amplitude has the form 
f d 2 qid 2 q'd 2 qL £J (_ a )j(-??) _^ 

A 2 _ 2+ 2r ~ y ( 27r) ' gr r ll l jgja $i(gl,-,e + i)^(g1,-,^ + i)- (64) 
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The impact factors $i ; 2 are given in terms of the integrals over the Sudakov parameters a\ 
2q[pA/s, [3[ = 2q[pB/s from the reggeon-particle scattering amplitudes f\ : 2 



r—1 r i r—1 



M JL 2" - fi Jl 2" 



The tree expressions for the amplitudes fi j2 appearing in the planar diagrams in QCD are 
given below 

fi = 1 1 1 1 

h (pa - q [? (pa -ho- q [? - ( PA - E[ =1 ^ - E[=o ^) 2 (pa - ELi 4 - E[=o ^) 2 ' 



ik = ^_L i i i 

h (pb + gi) 2 (pb - w + ,1)2 - {PB + E[=1 ^ _ Y:t% h? (pb + ELi - e 2 :^ 2 ^) 2 ' 

where fco = pa', k 2r+ i = PB'- The additional factors Ii j2 contain effective reggeon vertices for 
the production and scattering of the gluons with the same helicity. They can be written in the 
multi-Regge kinematics as follows 

T _ t t g|+i(Q - Et=i 4 - Et=l k t) {t r 

1 t\ (q*-e£U*-e^*) tr' 



l'l+l(Q* + Et=l It* Et=l ^2r~t+l) J-J 



r 

a r , 



where Q = pa — PA', Q = Pb — Pb> and the Sudakov variables of the produced particles 
ai = 2kipA/s, 0i = 2kips/s are strongly ordered 

1 > \fii\ > |/3 2 |- > \/hk\, \ai\ < H < ...|a 2 fc| < 1. 

In the physical region, where the signs of the Sudakov parameters of momenta k\ alternate with 
the index I 

/3i, a 2r < ; fo, «2r-i > ; 3 , a 2r -2 < ; ... , 
which is equivalent to the following constraints on the invariants 

si < 0, s 2 < 0, s r < 0, s r+ i > 0, s r+2 < 0, s r+3 < 0, s 2r+1 < 0, s > , (66) 

the integrands in expressions for <J>i 2 contain poles over the variables a' t , j3[ above and below 
the integration contours L over. Therefore are non-zero and can be calculated by taking 
residues from the poles 

■"' " r+l) = S (q* - eU 9? - e2 w - eS # - e2 a?) ' (67) 

*M, = n g - EU g , _ grr - gji - _ gzi ^_ t+i) • (68) 

In the case of production of 2r gluons with the same helicity the amplitude in N = 4 
SUSY is proportional to the Born expression. In the leading logarithmic approximation for the 
r + 1-reggeon contribution to the s r +i-channel the proportionality factor has the form 

Illa -y—) QQ J fiw^w^fi^p — ^7^ $l$2 ' (69) 
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where we introduce the new notation pi for the reggeon momenta q[. The multi-reggeon Green 
function satisfies the equation [35J 



l-2e 



G(p,^; Sr+1 ) = irG(p,p'; Sr+1 ), G(ftjf;0) = ]J^—S 2 - 2e ( Pl -p{) . (70) 

Here the kernel K in LLA can be expressed in terms of the infraredly stable Hamiltonian H 
K = u(t)-^H, w (i)=aQ -In -£) , t = -\q\\ (71) 

H = J P ^ P ^ 2 + ± Hl , l+1 , (72) 

where the pair Hamiltonian is 

Hi I+1 = In |p, | 2 + In |p m | 2 + pj" +1 In 1 2 — — + p[ p i+1 In \pi i+1 1 2 — — . (73) 

8 Integrable open Heisenberg spin chain 

The Hamiltonian for the gluon composite state in the ajoint representation has the property of 
the holomorphic separability [35J 



H = h + h*, fr = ln^±i + J> M+lj (74) 



Pi Pr+1 

■ -I- 

1=1 

where 



hi i+i = Inpi + lnpi + i + pi lnp/ ; + i — + pi+i lnp M+ i — !— . (75) 

Pi Pi+i 



Using the duality transformations (cf. [8]) 



d 

Pi = z 0,l , Pr = Zr-l,r , <? = Z0,n , Pr,r+1 = = id r , (76) 



the holomorphic hamiltonian can be rewritten as follows 

ra-1 

ZO.l Zn-.i. 



h = In — '—^ — — + ^2 K,r+i , (77) 

where 



Z 0,n r= i 



h r , r+ i = 2 ln(d r ) + — — + — — + +2 7 . (78) 

Here and later we neglect the pure imaginary contribution 21n(i) because it is cancelled in the 
total hamiltonian H. 

One can verify, that in the new variables h is invariant under the Mobius transformations 

azk + b 

z k -» — . 79 

cz k + d 

Therefore we can put 

z = , z n = oo, (80) 
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Further, by regrouping the terms one can write the holomorphic hamiltonian for n-reggeon 
interactions in the adjoint representation in other form [35j 



n-2 



h = -2\n z , n + ln(4i^) + ln (4-l,rA-i) + 2 7 + h' rr+l , (81) 



r=l 



where 



K,r+1 = ^(^r+l^r) + M^r.r+lA+l) ~ 2 In Z r , r+1 + 2 7 

= ln(d r ) + ln(9 r+ i) + -w- In z r r+ i <9 r + ^— In 2^+1 <9 r+1 + 27 . (82) 

C r GV+1 

The pair hamiltonian /i^, r+1 coincides in fact with the expression (j!4j) in the coordinate repre- 
sentation acting on the wave function with non-amputated propagators. 

The remarkable property of h is its commutation with the matrix element D(u) of the 
monodromy matrix (j!9)) introduced above for the description of the integrability of the BKP 
equations in the multi-color QCD 



[D(u),h]=0. (83) 
Therefore if we write D(u) as a polynomial in u 



71—1 



D(u) = J2^~ 1 ~ k 1k, (84) 



fc=0 



then the differential operators 



71-1 



1 , q[ = —i z r ®r > (85) 



r=l 

fe-l 



Qk 



Zri Y[zr s , rs+1 Y[idr t (86) 
0<ri<r2<— <rk<n s=l t=l 

are independent integrals of motion with the properties 

[q' k ,h] = [q' k ,q' t }=0. (87) 

It turns out, that h coincides with the local hamiltonian of the open integrable Heisenberg 
model in which spins are generators of the Mobius group. 

To solve this model one can use the algebraic Bethe anzatz. In this case it is convenient to 
go to the transposed space, where there exists the pseudo-vacuum state ^0 

71-1 

*o = II z r 2 > (88) 

r=l 

satisfying the equation 

C*(n)^ = 0. (89) 

Here C (u) is the transposed matrix element C{t) of the monodromy matrix (|19p . The eigenval- 
ues of the hamiltonian and the integral of motion D(u) are constructed by applying the product 
of its matrix elements £?*(«) to the pseudovacuum state 

k 

= n B *K)* . (90) 

r=l 
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For such eigenfunctions the spectral parameters u r should obey so-called Bethe equations. 
Instead one can introduce the Baxter function which is the generating function of the Bethe 
roots 

oo 

Q(u) = ~[[(u-u k ). (91) 

fc=l 

Generally the number of the roots Uf~ is infinite. The Baxter function satisfies the Baxter 
equation which is reduced to the simple recurrent relation for our open spin chain 

k{u)Q{u) = {u + i) n - l Q{u + i), (92) 

where A(u) is an eigenvalue of the integral of motion D(u) and can be written in terms of its 
roots 

71-1 

D{u) ^ ai ,a 2 ,.. . >an -i = A(«) *oi,aa,...,an-l » A (u) = ]J (« " *Or) ■ (93) 

r=l 

As a result, the solution of the Baxter equation can be found in the form [35] 

r=l ' 

up to a possible factor being a periodic function of — iu. 

The Regge trajectory of the composite state of n — 1 gluons has the additivity property 

w n (t) = w(t)-^, £ = e + ?, (95) 

n— 1 n— 1 

(: = ^2^{ar),e = ^2e(a r ), (96) 

r=l r=l 

where 

ft 

e(o) = ^(a) + V(l - a) - 2^(1) , a r = ^ r + y . (97) 

9 Three gluon composite state 

The wave funcion of the three gluon composite state in the adjoint representation can be 
constructed as a bilinear combination of eigenfunctions of the integrals of motion D(u) and 
D*(u) having the property of single- valuedness in the coordinate space |35j 

/j2 / 1 \ ai / * l\ ai 

11 !,-««( y*)-* (V—L) ,x = ^. (98) 

\y\ 2 \y-xj \y*-x*j z\ 

One can perform its Fourie transformation to the momentum space 

**(pi,A) = (pi+P2r ai - a2 (pi +pir~ ai ~~ a2 m ,y = -, (99) 

pi 



where 



m = J d 2 t + 1) " + 1) " (i - 1)-- 1 (i - ef' 1 ■ (ioo) 

This function can be presented in terms of its Mellin transformation 

// \ —iu / \ — iu 
d 2 u4>(u,u) (JiJ (j|J , (101) 
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where 



and 



iu = iv u + — ^ , -iu = iv u , / d 2 u = du u ^ . (102) 



7V u =-oo 



_ 7r 2 r(l + a l )T{a 2 ) T(iu)T(l+iu) T{-iu - ai)T{-iu - a 2 ) , , 

<^u, u) - r( _ a ^ r(i _ ^ r( _ r(i _ + _ + ^ )r(i + ~ + ■ i Udj 

Really the last form of \&* corresponds to the Baxter-Sklyanin representation [TT], because the 
function <fi is a product of the pseudovacuum state and the Baxter function [35] 

4>(u,u) = uuQ{u,u) , (104) 

where _ 

r(iu)r(m) r(-m-ai)r(-w-a 2 ) , , 

^ l ' ; ~ r(i - in) r(i - iu) r(i + »« + oi)r(i + »5 + a 2 ) ' { ' 



10 Discussion of obtained results 

It was demonstated, that Pomeron in QCD is a composite state of reggeized gluons. The 
BFKL dynamics is integrable in LLA. In the next-to-leading approximation in N = 4 SUSY 
the equation for the Pomeron wave function has remarkable properties including the analyticity 
in the conformal spin n and the maximal transcendentality. In this model the BFKL Pomeron 
coincides with the reggeized graviton. The BDS ansatz for scattering amplitudes in iV = 4 
SUSY does not agree with the BFKL approach in the multi-Regge kinematics. The reason 
for this drawback is the absence of the Mandelstam cuts. The BFKL-like equation for the 
composite state of two reggeized gluons with adjoint quantum numbers is explicitely solved. It 
is shown, that the equation for the composite state of an arbitrary number of reggeized gluons 
in the adjoint representation is equivalent to the Schrodinger equation for an integrable open 
Heisenberg spin chain. The wave function for three gluon composite state is constructed in the 
Baxter-Sklyanin representation. 

In the conclusion I thank L.D. Faddeev, J. Bartels and A. Sabio Vera for helpful discussions. 
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